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We derive kinetic equations describing injection and transport of spin polarized carriers in organic 
semiconductors with hopping conductivity via an impurity level. The model predicts a strongly volt¬ 
age dependent magnetoresistance, dehned as resistance variation between devices with parallel and 
antiparallel electrode magnetizations (spin valve effect). The voltage dependence of the magnetore¬ 
sistance splits into three distinct regimes. The hrst regime matches well known inorganic spintronic 
regimes, corresponding to barrier controlled spin injection or the well known conductivity mismatch 
case. The second regime at intermediate voltages corresponds to strongly suppressed magnetore¬ 
sistance. The third regime develops at higher voltages and accounts for a novel paradigm. It is 
promoted by the strong non-linearity in the charge transport which strength is characterized by 
the dimensionless parameter ell/ksT. This nonlinearity, depending on device conditions, can lead 
to both significant enhancement or to exponential suppression of the spin valve effect in organic 
devices. We believe that these predictions are valid beyond the case of organic semiconductors and 
should be considered for any material characterized by strongly non-linear charge transport. 

PACS numbers: 72.20.Ee 72.25.Dc 72.25.Hg 72.80.Le 


I. INTRODUCTION 

Magnetoresistance in organic spintronic devices is an 
easy effect to detect [I|. Vertical multilayered devices 
combining manganite (Lao.rSro.aMnOa) and cobalt ex¬ 
ternal magnetic electrodes with an Alqa transport inter¬ 
layer have perhaps never failed (no negative reports avail¬ 
ably to exhibit a clearly measurable magnetoresistance 
[2-l6j. This magnetoresistance, often called in literature 
giant magnetoresistance (GMR) or spin valve effect (SV), 
is generally measured as a difference between the device 
resistances for antiparallel and parallel electrode magne¬ 
tizations, where the latter is expected to facilitate the 
transport of spin polarized currents. Manganite based 
devices commonly show an inverse spin valve effect fea¬ 
turing lower resistance for antiparallel orientation. This 
property was explained recently by the spin filtering ef¬ 
fects at the organic-ferromagnetic interface Substi¬ 
tuting the manganite with a 3d metal or other magnetic 
electrode leads to a more controversial picture merging 
successes and failures. Generally, this second case has 
well established achievements in the tunneling regime, 
that is for a few nm thick layer of various organic semi¬ 
conductors Eini- On the other hand, the application 
of two metallic magnetic electrodes in combination to 
thick organic layers, expected to transport carriers via 
molecular electronic states, has lead to a number of re¬ 
ports claiming no magnetoresistance(see for example fl^ 
and high thicknesses regimes i n fiol ini) and a number 
of positive communications [ij, Il4| . 

An important fact to underline is the voltage depen¬ 
dence of the GMR in organic based devices. The mag¬ 
netoresistance is generally detectable at relatively small 
voltages (< 0.1 V) which are much lower than interfacial 
electronic barriers for the injection into two transport 


states of OSG, HOMO and LUMO. Indeed, typical values 
of the barriers measured by various spectroscopic tech¬ 
niques are of the order of 1 eV This clearly implies 
the need to consider eventual intragap states or band, 
as it is discussed in [l^, [l^ . It is also known from ex¬ 
periments that impurity states can be created in OSC in 
uncontrollable way during sample preparation [ 13 . 

Such easily achievable spintronic effects in electrically 
operated organic devices looked strange especially when 
compared to inorganic counterpart, where the hght with 
materials and concepts impeded similar achievements for 
years of tough and extremely qualified research. For in¬ 
organic devices these difficulties were convincingly ex¬ 
plained through the conductivity mismatch problem [3- 
[20|. A fundamental obstacle for the spin injection from a 
ferromagnetic metal to a semiconductor was formulated 
by Schmidt et. al.[T^. It was shown that in the case of 
diffusive transport the large difference in conductivities 
of ferromagnetic metal and semiconductor limits the spin 
injection coefficient 7 oc CTs/uf, where as,crF are conduc¬ 
tivities of semiconductor and ferromagnetic metal respec¬ 
tively. This limitation can be overcome by introducing 
a tunneling barrier at the ferromagnetic-semiconductor 
interface pol . [^ . leading to an essential improvement 
of the spin injection into a semiconductor. In order to 
be effective the tunnel resistance should be of the or¬ 
der of the effective resistance defined as Lp jap, Lsl^rs, 
where Lp,Ls are the spin diffusion lengths in ferromag¬ 
netic metal and semiconductor respectively. 

The question whether OSC are free of conductivity 
mismatch or the MR in OSC may not be caused by 
real spin injection into the organic electronic states has 
raised. The latter doubt was strongly enhanced by the 
absence of Hanle effect in all the tested organic spin¬ 
tronic devices [H, [U. Indeed in a device based on the 


2 


transport of spin polarized carriers, the application of a 
magnetic field perpendicular to the carrier’s spin causes 
the latter to process with a frequency 28 GHz per Tesla, 
regardless of whether the transporting medium is inor¬ 
ganic or organic. This precession causes a misalignment 
between the relative orientation of the injected spin po¬ 
larization and the electrodes’ magnetization, resulting in 
a measurable resistance change. Noteworthy, in inorganic 
spintronics Hanle effect has acquired the role of decisive 
proof for the demonstration of a spin polarized electrical 
injection in the investigated medium 

While interesting and stimulating ideas have already 
been advanced to explain the absence of the Hanle ob¬ 
servation we would like to underline an important 
paradox related to organic spintronic devices. Although 
magnetoresistance represents the parameter whose vari¬ 
ation is used to indicate the strength of the Hanle effect, 
there exist not yet a clear understanding of the nature of 
magnetoresistance in such devices. 

At present there were few atteimts to model spin 
transport in organic semiconductors [A . The first 

successful approaches to describe magneto-resistance in 
OSC were based on the extension of the models describ¬ 
ing inorganic devices [1^ through an accurate inclu¬ 
sion of the mechanisms of charge injection in OSC 23. 
An important step forward was proposed in Refs. [2^, 

The authors have considered a model of multistep tun¬ 
neling through single molecular level in OSC and the fi¬ 
nite occupation of the level during tunneling process. It 
promoted thus for the hrst time elements of nonlinear¬ 
ity for the spin polarized carrier transport inside OSC. 
Nevertheless, non linear effects were considered only for 
a small number of intermediate states, while for the bulk 
hopping case the resistance of organic was described as 
a constant Refs. [2^. 

Here we make a step further and consider charge trans¬ 
port nonlineartiy in an organic semiconductor confined 
between two spin polarized electrodes. We limit our 
model to the conductivity via a single intragap electronic 
level. In spite of its obvious simplification (compared to 
strongly non-homogeneous situation in OSC) the model 
clearly confirms two well known spintronic regimes de¬ 
scribed above, and namely the conductivity mismatch 
M and the tunnel barrier regime [iSlil, as expected 
from both experimental results and theoretical predic¬ 
tions. Surprisingly we discover a third regime, explicitly 
characteristic for organic semiconductors, where non lin¬ 
ear effects induce a measurable magneto-resistance. In¬ 
deed, for strongly voltage dependent hopping transport 
through the molecular level, the conductivity becomes 
strongly dependent on the injected spin. We demonstrate 
that in high voltage limit ell/ksT >> 1 the resistance of 
the organic semiconductor is governed by the contacts. 
Here e is the elementary charge, U is the voltage, fcs is 
the Boltzmann constant, and T is the temperature. In 
that case boundary conditions define the resistivity of 
organic semiconductor and therefore the resistivity be¬ 
comes spin dependent. It means that contrary to the case 


of ordinary semiconductor with diffusive transport the 
conductivity mismatch arguments are applicable for or¬ 
ganic semiconductors only for some limited voltage range, 
while for higher voltages strong nonlinear effects in or¬ 
ganic semiconductor enable the detection of the magne¬ 
toresistance caused by spin polarized injection. 

Although we have formulated the model for the con¬ 
ductivity via the impurity level, it is valid for any en¬ 
ergy level with hopping conductivity when the space- 
charge may be neglected. As a result the model is valid 
if the conductivity occurs over the HOMO and LUMO 
provided that concentration of the injected carriers is 
small. Experimental situation as far as transport in SV 
devices is concerned remains controversial. The trans¬ 
port due to electron and hole injection into the HOMO 
and LUMO is commonly assumed in literaturefl^ [Slj. 
On the other hand the I-V characteristics of the device 
in that case should be highly nonlinear and sensitive to 
the temperature (T^. This was not reported experimen¬ 
tally. These arguments lead to the suggestion that the 
conductivity of the SV devices is due to existence of some 
impurity levels when the applied voltages are much less 
than Fermi level-LUMO (and Fremi level - HOMO) dis¬ 
tance [l5l| . The origin of the impurity levels may be dif¬ 
ferent (for example metal atoms, oxygen molecules [l^ 
or X-rays induced traps during electron beam deposition 
of metallic electrodes [l^) . Contrary to Ref.flq where 
it was suggested that impurity levels have wide energy 
distribution between the HOMO and LUMO we assume 
that the impurity level has narrow energy distribution. 
We show in our work that for such levels the non-linearity 
can both suppress but also enhance spin valve effect in 
some cases. 


II. MAIN EQUATIONS 

In order to describe nonlinear effects we assume that 
the transport in an organic semiconductor in the low volt¬ 
age limit {ell << A, A is the LUMO-HOMO splitting) is 
determined by the polaron hopping over impurity levels 
similar to Ref. [T7I| . In order to have reasonable conduc¬ 
tivity the level should have energy close to the Fermi 
energy. We assume that the density of this levels n is rel¬ 
atively small (< 10^®cm“^)and therefore we can neglect 
the electric field which appears due to finite occupancy of 
the impurity levels by polarons. The electric field inside 
the semiconductor is equal to the external field in that 
case. We also assume that spatial energy fluctuations of 
these impurity levels are small Ae << fc^T. 

To describe the polaron hopping transport we apply 
the stationary kinetic equation for the distribution func¬ 
tion faim) in the site representation derived in Ref.js^: 

0 = y^^(fa{m'){l - faim)))Wm',m 

m' 

- fa{m){l-faim'))Wm,m'^ ( 1 ) 
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where 


Wm',m = exp {f}eE{m' — m))w(\m' — m\) (2) 

and w{\n\) = = ^/^bT, Ea 

is activation energy of polaron, h is the Planck constant, 
e is the charge of electron, and J{\n\) overlap integral 
between two sites. This equation are derived under con¬ 
dition that the spin relaxation is absent and the Hubbard 
energy is neglected. We consider only the case of weak 
field l3eEa << 1 (a is the hopping distance), therefore we 
can keep only lowest order in the expansion of the hop¬ 
ping probability in electric field E: Wm',m = w{\m' — 
m|)(l -|- PeE(rn' — m)). Assuming that distribution func¬ 
tion faim) is smooth function of m, we may expand 

fa{rn') = fa{m) + {Tn'-m) fa (m) + ^f^{m) 

and consider m as a continuous variable x. In that case 
Eq.(IT|) reads: 

D^-^+PeED-^{fa{x)il-fa{x)))=0 (3) 

Here the spin independent diffusion coefficient is defined 
as D = J2n w{n)n^a^/6. f3eDfa{l — fa) is the local con¬ 
ductivity of polarons with the spin a. The first integrals 
of these equations ( [S]) define the currents of spin ” up” 
and spin ’’down” polarons: 

= D^^+PeEDf^ix)il-Mx)) 
en ax 

— = + /3eEDf^{x){l - f^{x)) (4) 

en ax 


Differently from Ref. (2^ , in our case it is important to 
consider not only the change of the electro-chemical po¬ 
tential but also the contribution to it from the 

contact voltage. Note that the organic side of the contact 
cannot have surface charge, therefore the electric field in 
the contact is the same as in the bulk of organic layer. It 
means that we can neglect the contact voltage provided 
that contact is much shorter then the thickness of the 
organic layer. Therefore the change of the electrochem¬ 
ical potential is equal to the change of chemical 

potential. It means that the boundary conditions for the 
distribution function at the interface are defined as fol¬ 
lows: 


/Td(0) 

ftJL) 


I 

exp (/3(fo - J-t,;)) -I-1 ’ 

1 

exp {P{£o + + 1 


( 5 ) 


Here Sq is the position of the energy level in organic 
semiconductor with respect to the Fermi level. As a 
results Eq.(|4]) define distribution function f^,^{x) inside 
the organic layer. The substitution of this solution to 
the boundary conditions Eq.(IS]) allows us to obtain the 
equation for the current 

In order to discuss the current-voltage characteristics 
and the spin-valve effect we rewrite equations dH) in di¬ 
mensionless form: 

^^^yi^ + /Td(2/)(l-/td(2/)) =Jt.4./C (6) 


Here n is density of conducting levels. These equations 
are similar to ordinary diffusion equation describing spin- 
polarized currents in ordinary semiconductors [331 l34| . 
The main difference is in the dependence of the drift 
current on the distribution function. If the distribution 
function is small the current is defined by / << I. If 
the level is filled by electrons the current is determined 
by the concentration of holes (I — /). The later term is 
missing in the ordinary description of the spin transport 
in semiconductors described by ordinary diffusion equa¬ 
tion. These equations do not include the spin relaxation 
term (its effect will be discussed below). This holds if 
the thickness of the organic layer L < ^/Drf where Ts is 
the spin relaxation time in organic semiconductor. In the 
case of organic semiconductor ~ 10“^ — is 

large and ^/Dts is of the order of few hundred nanome¬ 
ters. 

Eqs. (HD should be supplied by the boundary condi¬ 
tions, which describe the microscopic structure of the 
contacts. Here we consider an oversimplified model (sim¬ 
ilar to that proposed in Ref. [2^) that assumes the ab¬ 
sence of the spin relaxation in the contacts. The change 
of electro-chemical potential in the contact is approx¬ 
imate by the Ohm’s law = JaRi,r-,a- Here index 

r, I stands for right and left interface and <t =t,4. stands 
for ”up” and ’’down” spin polarization. Ri^r,a are the 
contact resistivities for up and down spin. 


where dimensionless current defined as = 

J^^^L/enD, and dimensionless voltage f = (3eEL, y = 
ax, a = Q/L General solution of Eq.® is: 


h-xiy) = \- 2 tanh(K(;/ - yof.t)) (7) 

Here = ^(1 -b k), / 2 t,; = ^(1 - k) and k = 

(1 — 4j-|-,^/C)^^^- At the fixed voltage f currents and 
parameter are defined from the boundary condi¬ 

tions: 


/Td(0) = /ifd 

ftxiO = futX 


1 

exp(eo - Dt.tJtd)) + 1’ 

I 

exp (eo -b Mditd)) + 1 


( 8 ) 


where we introduce the dimensionless position of the 
energy level in organic semiconductor I3£q = Cq and 
the dimensionless contact resistances defined as = 

Rr,it,i/^o and the resistance of organic layer is Rq = 
L/fie^nD. Unfortunately the accurate solutions of these 
equations is possible only by numerical techniques. Some 
accurate analytical solutions are possible only in some 
limiting cases. The spin valve magneto-resistance is de¬ 
fined as follows: 


MR = 


jp - JA 
JA + jp 


(9) 
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Therefore in the case when total resistance is governed 
by the resistance of the organic layer we have MR = 

<C 1 in agreement with Refs.[ll, [l^. For 
large enough contact resistances Eq. (Ha gives a measur¬ 
able ma gnet oresistance in agreement with the results of 

Refs-liOnil. 

The spin-valve magneto-resistance decreases with ap¬ 
plied voltage, because effective resistivity increases with 
U. In the limit of intermediate voltage ell ^ ksT (( ^ 
1) the nonlinear effects become important. The right- 
hand side of Eq. cni) is exponentially small exp(—C) <C 1 
and may be neglected. In that case Eq. m has two 
roots leading to the transcendental equations for the cur¬ 
rent: 


FIG. 1: Spin-valve magneto-resistance MR as a function of 
dimensionless voltage eo = 2, r-|- = 1 r;. = 5. Organic layer 
resistance ro = 9.52. The low volta ge r egion (I) corresponds 
to the standard theory Refs.fl^. [^. l22l|. In the intermediate 
voltage region (II) the spin-valve is exponentially reduced. 
Region III corresponds to recovered spin-valve effect in the 
strong voltage region. 


where jA,p is the current of the device with parallel and 
antiparallel polarization of the ferromagnetic electrodes 
with the fixed voltage. For practical purposes it is conve¬ 
nient to exclude exp(—j/ot.t) from Eq.(|S]). As result we 
obtain a single transcendental equation for the current j: 

(/ltd “ - /zT.i) = 

e“^(/it.i - /fltd)(/it,4. - (10) 

Here C, = kC,. 


III. MAGNETORESISTANCE 


The results of numerical solution of Eqs. (nni) for cur¬ 
rents and spin valve magneto-resistance Eq.dH]) are plot¬ 
ted in Fig. 1. As it follows from this figure the spin-valve 
magneto-resistance has three characteristic regions. At 
small voltages eU « ksT we expand in Eq. (nni) /l, Rt,i 
up to the first order in 1- After straight¬ 

forward calculations we find: 


Jt.i = 


c 

nt.i + Md + ’’o 


( 11 ) 


where the resistance of organic layer is defined as = 
4cosh^ (eo/2)- Therefore in the limit of low field spin- 
valve magneto-resistance MR is governed by the conduc¬ 
tivity mismatch condition and by the contact resistances. 
Indeed in this limit we are able to reproduce well known 
results for ordinary semiconductors [isl [l^. In the 
symmetric case = ’’td spin-valve magneto¬ 

resistance is given by the formula: 


MR = 


_(?^t - _ 

(j't + + ^o)^ + (2r^ + ro){2r^ + r^)' 


( 12 ) 


/ltd ~ /i-td ^ j/td 


_c_ 

exp (-eo + I’/td/itd) + 1 


/2td = /flfd ^ /fltd 


_C_ 

exp (eo + rrt.tJiJtd) + 1 


They have two solutions jTtd JRtA respectively. 
The correct solution should satisfy the condition > 
/itd > /fltd — /ltd- These inequalities are satisfied 
if Jtd = ™*ii(ji?tdi/itd)- This formula defines current 
voltage characteristic in the large current limit. 

For the intermediate field case, when ksT <C eC/ fo 
(1 C Co), the correct root is determined by the 
sign of eo. In this case the conductivity is determined 
by the contact with positve energetic mismatch and the 
device shows no spin-valve effect. The actual value of the 
spin-valve magneto-resistance MR is exponentially small 
function of the applied voltage exp(—C). 

However at higher voltages when ^ the cor¬ 

rect root is determined not by the position of the organic 
level eo, but by the largest of the contact resistances. In 
the antiparallel (AP) configuration the root is always de¬ 
fined by the resistance r^. Here we again assume that 
r-j- > rj,. On the other hand in the parallel (P) con¬ 
figuration for one of spin directions the boundary with 
the resistance is absent and the root is determined by 
r^. It leads to reappearance of the spin-valve magneto¬ 
resistance at large voltage. In order to clarify these pro¬ 
cesses we plot in Fig.2 and Fig. 3 the spatial dependence 
of the level occupancy and the chemical potential in the 
organic level for the parallel and antiparallel configura¬ 
tions. As it is clearly seen, the resistivity of the device 
is determined by small density of holes in the left part of 
the organic layer and by the small density of electrons in 
the right part of the layer. Therefore the resistances of 
the antiparallel configuration is equal to the resistance of 
the parallel configuration P /• The parallel configuration 
P / has lower resistance leading to a large magnetoresis¬ 
tance. 

With the further increase of the voltage the equation 
for current may be solved with the logarithmic accu¬ 
racy ^ In (C). In this limit the spin-valve 













5 



FIG. 2: Spatial dependence of the level occupancy for spin up 
t and spin down electrons in the device with parallel ”P” 
and antiparallel ”A” configurations of contacts. 



FIG. 4: Spin-valve magneto-resistance MR as a function of 
dimensionless voltage ^ for = 1 r^, = ^ and different posi¬ 
tions of conduction level cq. 



FIG. 3: Spatial dependence of the chemical potential for spin 
up t and spin down in the device with parallel ”P” and 
antiparallel ”A” configurations of contacts. 


magneto-resistance increases, logarithmically approach¬ 
ing the limiting value: 


MR = 


(r^ -h Sr;) 


(13) 


which is larger then weak field value (fT^ This formula 
does not contain any information about the resistance 
of the organic layer and therefore represents the perfect 
device. 

We have described above three limiting cases for the 
spin valve effect. It is important to underline that the 
position of the conducting level cq with respect to the 
chemical potential of ferromagnetic contacts is a key pa¬ 
rameter for establishing a magnetoresistance (spin valve 
effect). Figure 4 shows the dependence of the magne¬ 
toresistance on the difference between level energy and 
chemical potential of the ferromagnetic layer eg for con¬ 


stant resistances. When Cq < 1 the spin-valve effect 
monotonically grows with applied voltage reaching the 
asymptotic value Ea. dO)) (Fig.4). When ep = 1.5 the 
minimum of spin-valve effect at intermediate valuees of 
electric fields Q. At ep = 2 we clearly see the three re¬ 
gions, including that where the spin valve effect is fully 
suppressed. 


IV. SPIN RELAXATION 

Neglecting the spin relaxation (as it is done in Ref. p3 | ) 
it is possible to suppress conductivity missmatch only by 
introducing the contact resistances. When the contact 
resistances become larger than the resistance of the nor¬ 
mal layer, the device is not subject of the conductivity 
mismatch and high SV effect can be achieved. 

This is not the case if even a weak spin relaxation is 
included. It was shown in [23,[2l| that large contact resis¬ 
tances lead to strong spin accumulation and depolariza¬ 
tion near the boundary. So the contact resistances should 
be selected in a narrow range to avoid both phenomena: 
spin depolarization and conductivity missmatch. 

Indeed our model differs from one discussed in 
and [ 2 ^. In our case the conductivity missmatch is sup¬ 
pressed not only by the contact resistance but also due 
to applied voltage and non-linearity. Also we have some 
limit for spin accumulation because in our case the num¬ 
ber of electronic states is finite and we can not have 
or larger than unity. 

Let us now include in our equations a spin relaxation 
time T that does not depend on applied voltage 

D^^^+l3eED^{Ux)il-U{x))) = 

dx^ dx ^ ' T 


After introducing spin currents in the dimensionless 
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form we obtain 


dy 


dy 



/i), 7 


Dr’ 




/t) = Jt/C 


dy 


+ /t(l 


fi) = JiK 


This set of equations is much more complex then equa¬ 
tions without spin relaxation and allows only numerical 
analysis. Let us also note that the perturbation theory 
for small 7 is not useful in our case because the correct 
perturbation parameter appears to be and quickly 
becomes large at hight voltages. Therefore we solve these 
equations numerically for the case eo = 0 . 



FIG. 5: The magnetoresistance for the device with different 
spin relaxation rates 7 . Line correspond to analytical solution 
of the differential equations with 7 = 0 . 


Fig. [ 5 ] shows the numerical solution for a device with 
r-|- = 5, = 1, ep = 0 and different relaxation rates 

7 . The results for 7 = 0 agree with analytical solution. 
It can be seen that relatively small spin relaxation 7 = 
0.1 does not alter magnetoresistance significantly even 
at high voltages C- The tendency of magnetoresistance 
to grow with applied voltage remains even for the large 
relaxation rates 7 > 1 . 


V. DISCUSSION 

In a fully metallic device the conductivity of the normal 
layer does not depend on the shift of chemical potential. 
The description of the device in terms of electrochemical 
potential does not require the knowledge of the separate 
contributions from electrical and chemical potentials. In 
our case, where the conductivity is governed by few func¬ 
tional levels near Fermi level, the situation is different. 


Small shifts of the chemical potential change drastically 
the conductivity of the intermediate layer. Indeed, as¬ 
suming the existence of the local chemical potential via 
the distribution function 

f =_i_ 

expiPiSo-ya{x))) + l 
we may rewrite Eqs.(|3]) in the standard form 


:k 

en 



where is electrochemical potential and cr{fj,) = 
j3eDf^{l — fa-) is chemical potential dependent conduc¬ 
tivity. Therefore our theory allows standard formulation 
in terms of nonlinear conductivity. 

The shift of the chemical potential always occur in the 
spin-valve device with injection, although in devices with 
direct ferromagnet-normal layer contact (without contact 
resistance) it has slightly different form [l^, [l^ . There¬ 
fore the non-linear effects in magnetoresistance should be 
considered in all devices where the shift of the chemical 
potential of the order of applied voltage can significantly 
effect the conductivity. 

Let us consider the case when the transport in the or¬ 
ganic layer is governed by the HOMO and LUMO. The 
voltage ~ O.IV is not sufficient to inject carriers into 
these states and all the carriers in the HOMO and LUMO 
are thermally activated. The number of such carriers is 
small and it allows us to neglect the space charge mak¬ 
ing our theory applicable to this case. Transport via 
the HOMO and LUMO corresponds to the large energy 
£0 ~ leU and leads to the exponential suppression of 
the spin-valve effect with applied voltage when eV > T 
(see also [ 2 ^). This exponential suppression was never 
observed in the experiments. It gives another argument 
(besides the ones discussed in [Tg) against the HOMO 
or LUMO transport in the organic spin-valves. 

These arguments are in line with the arguments pre¬ 
sented in Ref. [T^ invoking the presence of impurity states 
in the organic gap to describe the organic SV devices. 
The nature of these states however is far from being 
clear. It was reported [13 that similar defects may be 
created during sample preparation, leading to the sup¬ 
pression of the spin valve effect in sufficiently thick sam¬ 
ples. The concentration of such defects can be very low 
(sub-percent level), precluding their detection by spec¬ 
troscopic techniques [ 13 . The author of Ref. [l^ consid¬ 
ered impurity states with concentration ^ 10^®cm“^ and 
broad distribution of energies. Such a broad distribution 
yields cr(/i) Ri const and the discussed non-linearity is 
absent. It was also assumed [l^ that the organic layer is 
highly conductive and its resistance is negligible relative 
to contact resistances. Note that the conductivity of a 
system of localized states with a broad energy distribu¬ 
tion is governed by the small number of states with the 
energy close to the Fermi level [s^ . Therefore the actual 
concentration of impurities that contribute to conduc¬ 
tivity is ~ which can hardly lead to the high 

conductivity of the organic layer. 
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There are two possible solutions of this problem. The 
impurity states may be arranged in some sort of conduct¬ 
ing channels or their energy distribution may be more 
narrow than assumed in Ref.[l^. The second option 
leads to the strong non-linear phenomena related to cr(/r) 
dependance. In the present paper we discussed these phe¬ 
nomena in the extreme case when the width of energy 
distribution of impurity states is less than the tempera¬ 
ture. 


VI. CONCLUSION 

We derived a set of kinetic equations describing non 
linear effects related to injection and transport of spin 
polarized carriers in organic semiconductors with hop¬ 
ping conductivity. The model predicts a strongly voltage 
dependent magnetoresistance splitted into three distinct 
regimes. The first regime (low voltages) corresponds to 


the well known in inorganic spintronics conductivity mis¬ 
match limitation, underlining the correctness of the ap¬ 
plied approach. The second regime at intermediate volt¬ 
ages corresponds to fully suppressed magneto-resistance. 
Interestingly, a third regime develops at higher voltages 
and accounts for a novel and purely organic paradigm. 
It is promoted by strongly non-linear effects in organic 
semiconductor which strength is characterized by the di¬ 
mensionless parameter eU/ksT. This nonlinearity, de¬ 
pending on device conditions, can lead to both significant 
enhancement or to exponential suppression of the spin- 
valve effect in organic devices. Consequently, while in in¬ 
organic devices the conductivity mismatch limitation can 
be lifted off only by inserting a spin dependent tunnelling 
resistance, organic devices feature an additional mecha¬ 
nism based on nonlinear charge transport behaviour. 
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